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^ ■ Abstract 
O ; 

\ Bose condensation of interacting bosons in a two-dimensional random potential 

^ \ is studied. The Gross-Pitaevskii equation is solved to determine the spatially-varying 

OO ' order parameter and the localization length as a function of the disorder, the interaction 

. strength, and the condensate density. A finite temperature of condensation is obtained 

\ thereby. The results are applied to determination of the superradiant decay of excitons 

in a GaAs quantum well. 

> ' 
C ■ 

; PACS number(s): 03.75.Fi, 71.23.An, 75.40.Cx 

^ • Bose condensation in a disordered medium was the topic of extensive research 

starting a decade ago |[l|, 0]. It was mainly concerned with the superfluidity of liquid 
^He in porous media or on substrates [^, ^, H] or with the superconductor-insulator 
^ ' transition in high critical temperature superconductors fs], 0, ^, ^. In most cases, the 

Hubbard Hamiltonian of interacting bosons on a lattice is employed. This model is 
treated via a mean-field theory via the Bogolyubov approximation [11 1, in the 



approximation of Gaussian fluctuations ||13|, by the Monte Carlo numerical methods 



|13, 14, 15], or by renormalization group methods [16, 17]. 

Bose condensation in other systems is often treated via the Gross-Pitaevskii equa- 
tion (nonlinear Schrodinger equation) for the condensate [p!^ ]. The advantages of this 
approach are that one immediately obtains the order parameter of a non-uniform con- 
densate and that the non-physical features related to the integer occupancy number 
of a lattice site do not appear. The connection between two approaches was made in 
Refs. i, 0]. 

Our interest in Bose condensation in a two-dimensional space is primarily moti- 
vated by the ongoing experiments aimed at observing quantum degeneracy of excitons 
in quantum wells [19[. Since in two or smaller dimensions the critical temperature 
of condensation would be zero in a uniform infinite space, only localization by ge- 
ometrical size, or, in our case, localization by disorder can cause condensation at a 
finite temperature. This phase, consisting of a number of separate and presumably 
phase-uncorrelated "lakes" of Bose condensate is more properly termed "Bose glass" 

i- 
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For excitons, the order parameter is proportional to the polarization of the medium. 
Therefore Bose condensation has a profound influence on the optical properties of 
excitons, especially on the superradiant decay. The smallest of the exciton mean free 
path and the localization length gives the coherence area of an exciton. It was shown 



in Ref. |20, 21, 22(1 that the coherence area is related to the linewidth and the radiative 



lifetime of excitons. Bjork et al. discussed it in terms of superradiance of excitons 



1 23]. A sharp increase of the exciton oscillator strength was suggested as a criterion 



of Bose condensation of excitons 24 1. The oscillator strength reaches saturation 
as the coherence length approaches the wavelength. Other approaches (Gerrit Bauer, 
private communication) predict a divergence of the coherence area as the temperature 
decreases down to some critical temperature, and a constant coherence area below this 
temperature. In view of this controversy, it is important to give an intuitive treatment 
of the oscillator strength and the superradiant decay. 

In this paper we numerically solve the Gross-Pitaevskii equation for the order pa- 
rameter of the condensate in a disordered two-dimensional medium modeled by a set 
of potential wells with their depths depending on the widths. [The authors understand 
the limited nature of such an assumption, but stress that such a potential might be 
applicable to quantum dots defined by the variation of a quantum- well width] . In this 
way we determine the localization size as a function of the disorder strength, the in- 
teraction strength, and the condensate density. This relation along with the integral 
of the Bose distribution for a localized condensate allow us to obtain the temperature 
of condensation for a given density of bosons. The localization length determines the 
oscillator strength and the increase of the rate of superradiant decay. 

We describe the condensate of interacting bosons via the stationary Gross-Pitaevskii 
equation for the long-range order parameter ip 

-^V^iJ + V{r)i; + Um^ =f^^, (1) 

where M is the mass of a boson, U is the interaction strength, fi is the chemical 
potential, and V is the external potential. This model implies that we take the lifetime 
of the boson particles to be infinite and neglect scattering. 

This would be a good approximation for indirect excitons (having the electron and 
the hole in different quantum wells) in a direct-bandgap semiconductors. Equation 
(|l|) also implies that the exciton-exciton interaction can be modeled by a short-range 
s-wave scattering. Even though this assumption is at best questionable, we make it for 
the sake of simplicity. 

To incorporate a random potential, caused, e.g., by impurities or fluctuations of 
the quantum well thickness, we adopt the following model. The whole two-dimensional 
space is divided into a set of "lakes" of condensate of the characteristic size Lc- The 
condensate in each of them is presumed to be trapped by a cylindrical potential well 
of a radius L to be determined: following Ref. [^] we choose the depth of this well to 
be inversely proportional to its radius 

y(r)=-|0(L-|r|), (2) 

where ^ is the parameter characterizing the strength of the disorder. The average 
density of bosons in the condensate is ric- The average number of bosons in each lake 
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is Nc = ricL^. We impose a condition to normalize the number of bosons in each 
condensate lake 

J iVpdr = n,Ll (3) 
and we identify Lc with the coherence length defined by the following relation 

J\^\'dr = nlLl (4) 

The integrals have infinite limits. 

In order to reveal the self-similarity of the solutions, let us introduce a constant 
length scale Lq and correspondingly express all densities in terms of On the 

other hand, the coordinates will be expressed in terms of the well width L and the 
wavefunction will be expressed in terms of 1/L. Then the Gross-Pitaevskii equation in 
dimensionless variables becomes 

- V^^ - ^oLd{l - \r\)i, + u^lV'P = fioL^^, (5) 

where the dimensionless parameters 



Co = — , (6) 



/^O - ^72 • V^) 



The parameter Ih is called "healing length" of the condensate, 

I' = = A fo) 

^ MUric uric' ^ ' 

By comparing the first and the third terms in (^), one can estimate that the order 
parameter can appreciably change only over the dimensionless length of the order of 
Ih- u is large, the healing length is less than the interparticle separation, and if u is 
small, it is larger that the interparticle separation. Under this change of variables (^) 
and (^) become 

" \^P\^dr = ncLl (10) 
\Ur = nlLlL\ (11) 
At any given strength of disorder ^, for convenience, we chose the length scale 

L - 

2M^ 



(12) 



so that ^0 = 1- Besides, we notice a self-similarity of the set of the equations (^)-(|T] 
under the transformation with an arbitrary positive a 



u' = u/a, n ^ n' = an, ip ^ ip' = ^/atp. (13) 
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Figure 1: Dependence of the chemical momentum of the condensate on the trapping well 
width at M = 1, n = 0.06. 

Therefore we numerically solve the set of equations for u = 1 and obtain the general 
solution by rescaling according to (|T3|). 

The solution of the dimensionless equations (|5|) - ( [111 ) determines the dependence 
of the chemical potential of the condensate on the trapping radius L. An example 
of it for M = 1 is shown in Fig. |l[ For the trapping radii outside the region on the 
plot no localized solutions (with /i < 0) were found. At u = 0, a linear Schrodinger 
equation, solutions exist for any radius L of the well. As the density of bosons, or, 
equivalently, the nonlinear coefficient u, increases, the interaction of bosons screens 
the random potential more strongly [^. As a result, the range of L, where localized 
solutions are possible, narrows until it totally disappears. The radius L is determined 
from the solution of (||) which minimizes the chemical potential fi. Note that extended 
condensate solutions are not a limit of localized solutions. The former have a positive 
chemical potential = uric- The latter must have negative n to ensure decrease ip ^ 
at |x| — > oo. 

Based on these solutions, we seek the dependence of the condensate size Lc on the 
average condensate density ric- The numerical result for n = 1 is shown in Fig. ^. With 
a good accuracy at all densities this dependence can be approximated by 

L, = a{ng - nf (14) 

with a = 5.4, [3 = —0.1317, and ng = 0.074. Thus below the critical density (propor- 
tional to the dimensionless constant n^), the condensate exists as a Bose glass with 
localized "lakes" of condensate, and, above this density, an extended Bose condensate 
is possible at zero temperature. This universal dependence will be used later while 
treating bosons with different nonlinear coefficient u and at non-zero temperatures. 
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Figure 2: Dependence of the localization length of the condensate on its density at u 

Now we apply our model to the case of excitons in GaAs quantum wells. Let us 
estimate the interaction strength U for excitons in the ideal two-dimensional case. It 
is proportional to the exciton-exciton interaction term, i.e., the Coulomb exchange 
integral in the exciton ground (Is) state, which, according to [^6| is equal to 



Q1 cJttS 

I = Y^lMk') - <Pis{k)Wu{mis{k)\'V{k - k') = 4^a2(l - -^)Eo, (15) 

kk' 

where (j)is{k) is the wavefunction of the Is exciton state and V{k) is the Coulomb 
potential in the momentum representation, qq is the Bohr radius and Eq is the binding 
energy of an exciton in two dimensions given in p^ . We have neglected screening of 
the Coulomb potential by other excitons, which is possible to do for small densities 
(interparticle separation much larger than oq). At higher densities (|l^) overestimates 
the strength of interaction u. 

The estimate of / for indirect excitons will be inserted here. 

Therefore we obtain that 

AMI 6.06M , , 

u = — K— ~ , (16) 

bearing in mind that M = me + uih and nir = niemh/M. For the case of GaAs, where 
rUe = 0.0665mo, rrih = 0.377mo and niQ is the mass of a free electron, we obtain u ~ 47, 
which means strong boson-boson interaction. It is the opposite of Bose condensation 
in atomic traps in which one encounters a weak interaction, 

u~- , (17) 

-'^trap 
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where as is the scattering length in s-wave colhsions, and Ltrap is the size of the trap. 

It is weh-known that the critical temperature of Bose condensation in two dimen- 
sions is zero for an infinite space. Therefore only a localized (Bose glass) rather than 
an extended condensate can exist at non-zero temperatures. Only a fraction (if any) 
of bosons will be in the condensate and be described by the Gross-Pitaevskii equation. 
1^ We need to add the relations describing the statistics of bosons outside the con- 




densate. Following 1 28 1 we obtain for a two-dimensional localized Bose gas with the 
average density n 

(18) 

nA,^, = log(^), (19) 

where ks is the Boltzmann constant, the de Broglie wavelength corresponding to the 
temperature T is 

and Acr is the wavelength corresponding to a critical temperature for the given lo- 
calization size Lc- Here we again assume the same scaling of lengths and densities. 
Equations (|T8|) and (|l9|) together with the rescaled ([l^) 

Lc = a{ng - uricY (21) 

allow us to determine the condensate density Uc and the localization length Lc for 
the given concentration and temperature. We see that the condensate density cannot 
exceed the critical value Ug/u. 

The correlation length allows us to obtain the collective dipole of the condensate 
and the rate of radiative decay. Exciton condensate is similar to a collective of excitons 
excited by a laser. The differences are: 1) that the direction and the phase of the dipole 
are not given by the laser, but undetermined until a measurement is made on the 
condensate; 2) the order parameter does not contain a spatial dependence associated 
with the wavevector of the exciting laser. With these modifications, we apply the 
results of Bjork et al. The number of modes available for excitons are 



o 2 

Ne = (22) 



On 



^0 

where the condensate radius is such that -Kai = Ll. In the case of a thin quantum well 



and large N^., according to [29|, the power of emission per unit angle is proportional to 



/((/), x)r (</>), where the usual dipole pattern 

^(((A, X) = cos^ X + sin^ X cos^ (23) 
and the factor of the collective emission is 

V kac sin cp J 



^We disregard the possibility of the Kosterliz-Thouless transition, see e.g. 
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Figure 3: Fraction of the bosons in the condensate (middle curve) the locahzation length 
in /im (lower curve) and the factor of enhancement of radiative decay rate (upper curve) as 
functions of temperature for u = 47 and Lq = 10~^m: a) n = 1.2 x 10^°cm~^. 

Here k = 2-k/X, and A is the wavelength in the material (228nm in GaAs), (j) is the angle 
between the direction of emission and the orthogonal to the well. The cooperativity 
parameter, which determines how well the radiators are localized, is 



/xc = /_ sin 0#(1 + cos" (/.)r(<^). (25) 







For small condensates {kac ^ 1) it tends to 1, and for large condensates it tends to 
3/{kac)^. The rate of superradiant emission of excitons 7 is increased compared to the 
bulk rate of free-electron-hole recombination 70 by the enhancement factor 

7 = /iciVeTO- (26) 

For small condensates the emission rate grows as to the same extent in all directions. 
As the size Oc of the condensate grows, the enhancement factor tends to a constant. 
However the angular distribution becomes strongly peaked at (j) = where it still grows 
as a^. 

The results of the numerical solution of (p!^)-(21) are shown in Fig. For illustra- 
tive purposes we restore the dimension of length by taking a specific disorder strength 
which corresponds to Lq = 10~^m. We see that the condensate (a non-vanishing share 
of the condensate density) appears below a well-defined temperature. If the total den- 
sity of bosons is below the critical density, the fraction of the condensate tends to 1 
as the temperature decreases, and the localization length tends to a small finite value. 
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Figure 4: Same as above: b) n = 0.8 x lO^^cm 
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Figure 5: Same as above: c) n = 0.4 x lO^^cm 
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see Fig. ^ If the total density is above the critical one, the density of the conden- 
sate tends to that critical value so that the condensate fraction is always less than 1. 
The localization length tends to infinity as the temperature decreases and its value 
strongly depends on the total density, see Fig. However this increase does not be- 
come noticeable until temperatures several times below the condensation temperature. 
The enhancement factor grows appreciably until it reaches the limiting value as per 
reasoning above. 

In conclusion, we have presented a theoretical description of Bose condensation 
in a random potential in two dimensions. The disorder is believed to be sufficiently 
strong so that the Kosterliz-Thouless transition the effects of the sample size can be 
disregarded. We argued then that only localized condensates (Bose glass) can exist 
at non-zero temperatures. A macroscopic occupation of the condensate appears in a 
sharp transition below some condensation temperature. The average density in the 
condensate always stays below a certain critical density determined by the disorder 
strength. As the temperature decreases, the localization length, and consequently the 
coherence area, grow indefinitely if the average density of bosons is higher than the 
critical density. The increase of the radiative decay rate reaches its limit when the 
correlation length reaches the photon wavelength in the medium. The increase in the 
rate of superradiant decay and, especially, in the power of emission orthogonal to the 
well, is a reliable indicator of Bose condensation. 

This work is supported by NSF, Packard Foundation, and the Center for Quantized 
Electronic Structures (QUEST), an NSF Science and Technology Center. We thank 
Leonid Butov for enlightening discussions and making available his experimental results 
before publication. 
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